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In  this  paper,  the  energy  management  problem  for  a  Fuel  Cell  Plug-in  Hybrid  Electric  Vehicle  is  inves¬ 
tigated  using  a  real-time  power  splitting  method  that  preserves  the  battery  pack  life.  The  reduction  of 
high-power  demand  during  the  vehicle  stop-and-go  operations  and  the  limitation  of  the  batteries’ 
minimum  energy  to  a  prescribed  level  are  the  two  aspects  considered  in  this  study.  Using  a  longitudinal 
motion  model  of  an  experimental  small  utility  truck,  a  multi-criteria  optimization  problem,  including  the 
hydrogen/electricity  cost  ratio  and  the  anticipation  of  the  high-power  demand,  is  formulated  and  solved. 
The  analysis  of  the  effects  of  the  fuel  cell  and  the  battery  efficiency  degradations  demonstrated  that  there 
exists  a  minimum  cost  ratio  that  prevents  the  Energy  Management  System  to  adequately  limit  the 
batteries’  depletion  with  the  contribution  of  the  Fuel  Cell  power.  The  overall  control  system  does  not 
require  a  prior  knowledge  of  the  trip  information.  Therefore,  the  real-time  implementation  of  this  system 
is  realistic.  Compared  with  a  rule-based  Energy  Management  System,  the  proposed  optimal  and  real¬ 
time  approach  provides  up  to  5%  of  hydrogen  saving  whilst  maintaining  the  Battery  Pack  energy  close 
to  the  minimum  prescribed  energy.  The  high-power  demand  reduction  and  the  depletion  limitation 
enable  the  Energy  Management  System  to  preserve  the  batteries’  life. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  use  of  fossil  energy  in  the  transportation  sector  is  one  of  the 
greenhouse  gas  emission  sources  [1].  So,  the  vehicle  powertrain 
electrification  and  hybridization  have  been  identified  as  a  way  to 
reduce  these  pollutants  and  to  increase  the  fuel  usage  efficiency  [2], 
Among  the  commonly  found  powertrain  electrification  approaches 
are  the  hybrid  electric  (HEV),  the  plug-in  hybrid  electric  (PHEV)  and 
the  all  range  electric  powertrains  which  are  based  solely  on  batteries 
(BEV)  [2],  Nowadays,  the  lack  of  wide  deployment  of  charging 
stations,  combined  with  their  limited  operating  range,  make  it 
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difficult  for  the  BEVs  to  achieve  the  same  travel  distance  as  the 
conventional  gasoline  vehicles.  The  PHEV  with  an  internal- 
combustion  engine  can  provide  an  extended  range.  However,  the 
use  of  fossil  energy  contributes  to  the  greenhouse  gas  emission  [3,4], 
Hydrogen  has  been  identified  as  a  clean  energetic  vector,  and 
therefore,  it  can  help  to  reduce  greenhouse  gas  emissions  when 
used  for  electrical  propulsion  [5],  Furthermore,  on  a  strictly  ener¬ 
getic  point  of  view,  the  hydrogen  specific  energy  (33.33  kWh  kg-1) 
is  by  greater  than  the  best  currently  available  battery  specific 
energy  (0.25  kWh  kg-1).  Hydrogen  can  be  used  through  a  fuel  cell 
or  through  a  genset  (engine-generator  powered  by  hydrogen)  [6], 
The  PHEV  serial  topology  allows  the  range  extender  power 
source  to  charge  the  batteries  [2],  Since  this  range  extender  power 
source  is  not  required  to  track  the  vehicle  power  demand  dynamics, 
it  can  be  set  to  operate  at  it  maximum  efficiency.  However,  the 
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Fig.  1.  Fuel  cell  hybrid  electric  vehicle:  the  arrows  with  the  solid  lines  represent  the 
energy  flows;  the  arrows  with  the  dashed  lines  represent  the  control  and  measured 
signals. 


battery  power  source  and  the  electric  traction  system  must  be 
powerful  enough  to  adequately  track  the  vehicle  user  power 
demand  [7,8], 

In  this  paper,  the  energy  management  and  the  power  splitting 
problems  are  thoroughly  analyzed  and  a  real-time  optimal  solution 
is  proposed  for  a  typical  fuel  cell  and  battery  plug-in  hybrid  electric 
vehicle  (FC-PHEV)  with  a  serial  topology.  In  addition,  we  focus 
mainly  on  the  Proton  Exchange  Membrane  Fuel  Cell  (PEMFC) 
because  of:  (i)  the  low  operating  temperature,  (ii)  the  possibility  for 
high  current  density  and  (iii)  the  good  behavior  under  dynamic 
operating  conditions  [9], 

It  is  well  documented  that  the  vehicle  performance  and  the  fuel 
savings  are  related  to  the  power  flow  control  between  each  energy 
source  and  the  powertrain  [10,11  ].  Among  reported  methods  for  the 
power  flow  control  of  FC-PHEV  are  the  rule-based  methods  [12], 
the  off-line  optimal  method  and  the  near-optimal  on-line  method 
[13],  The  rule-based  methods  use  a  set  of  predefined  rules  to  split 
the  power  demand  between  the  fuel  cell  and  the  batteries.  These 
methods  are  simple,  real-time,  implementable  and  robust  to  traffic 
perturbation.  Since  they  are  mostly  related  to  the  Energy 
Management  System  designer,  they  cannot  be  widely  used  [2],  In 
addition,  with  rule-based  methods,  there  is  no  formal  way  to  assess 
the  overall  optimality.  The  off-line  optimal  methods  aim  at 
providing  a  more  formal  way  to  design  the  FC-PHEV  energy 
management  and  power  splitting  systems  [14],  However,  their  non- 
causal  nature  limits  their  use  for  real-time  applications.  The  on-line 
optimal  methods  can  be  seen  as  the  real-time  version  of  the  off-line 
methods.  Usually,  a  local  cost  function  is  minimized  with  the  in-line 
method  whereas  a  global  cost  function  is  used  to  formulate  the  off¬ 
line  solution.  Some  of  the  reported  energy  management  methods 
have  used  the  off-line  solutions  in  order  to  adapt  the  in-line  solu¬ 
tion  when  the  vehicle  is  moving  [15,16], 

In  most  of  the  reported  works  on  the  fuel  cell  vehicles,  the  fuel 
cell  real  behavior  is  not  taken  into  account.  Nevertheless,  it  is  well 
documented  that  some  fuel  cell  dynamics  cannot  adequately  track 
the  vehicle  power  demand  [17,18],  Even  if  the  PEMFC  can  sustain 
the  power  demand  [19],  the  hybridization  with  a  battery  will 
improve  the  vehicle  efficiency  and  will  contribute  positively  to  the 
fuel  cell  life  span.  A  rapid  power  variation  imposed  on  the  fuel  cell 
may  result  in  the  oxygen  starvation  [20],  On  the  other  hands,  the 
Battery  Pack  is  sensitive  to  the  ambient  temperature  and  its  life 
span  is  related  to  the  depth  of  discharge,  the  number  of  charge/ 
discharge  cycles,  the  high-power  discharge  rate,  etc  [21  ]. 

This  paper  addresses  the  energy  management  and  the  power 
splitting  problem  with  two  main  contributions.  The  first  contri¬ 
bution  is  related  to  the  analysis  and  the  design  of  a  real-time  and 
efficient  Energy  Management  System  that  incorporates  the  effect 
of  the  different  sources’  efficiency  variation  over  the  time.  As  the 


PEMFC  Output  Power  (W) 


Fig.  2.  Experimental  Fuel  Cell  efficiency  and  power  production. 

FC-PHEV  is  operating,  various  factors  affect  the  power  source  effi¬ 
ciencies.  In  particular,  we  investigate  the  influence  of  the  battery 
power  source  efficiency  degradation  on  the  power  splitting 
method.  The  second  contribution  is  related  to  the  design  of  an 
optimal  power  splitting  algorithm  that  preserves  the  battery  life  by 
limiting  the  depletion  and  by  reducing  the  high-power  discharges 
originating  from  the  vehicle  stop-and-go  operations. 

The  rest  of  the  paper  is  organized  into  six  sections.  The  FC-PHEV 
control  architecture  and  the  vehicle  model  are  presented  in  Section 
2  and  3,  respectively.  Using  this  model,  the  energy  management 
problem  for  the  battery  life  preservation  is  formulated  in  Section  4 
and  solved  in  Section  5.  The  simulation  results  and  the  comparative 
study  are  discussed  in  Section  6,  whereas  the  conclusion  is  pre¬ 
sented  in  Section  7. 

2.  FC-PHEV  serial  control  architecture 

With  the  serial  topology  of  a  FC-PHEV  as  shown  in  Fig.  1,  the  fuel 
cell  can  operate  at  its  maximum  efficiency  when  recharging  the 
batteries  during  the  vehicle  operation.  We  assume  that  the  battery 
pack  is  adequately  sized  in  order  to  meet  the  power  demand  during 
a  trip.  Fig.  2  shows  an  example  of  a  real  small  scale  PEMFC  efficiency 
as  a  function  of  the  power  production.  In  this  case,  the  maximum 
efficiency  of  0.52  (see  Equation  (7))  is  obtained  with  an  output 
power  of  150  W. 

In  Fig.  1,  the  uni-directional  DC— DC  converter  is  designated  as 
the  fuel  cell  power  source,  whereas  the  battery  pack  and  the  bi¬ 
directional  DC— DC  converter  is  considered  as  the  battery  power 
source.  The  power  demand  Pu  from  the  vehicle  driver  is  given  to  the 
Energy  Management  System  (EMS)  as  a  sequence  of  set  points. 
Taking  into  account  the  available  energy  in  the  hydrogen  tank  EfC, 
the  available  energy  in  the  battery  pack  £b  and  the  different  power 
source  efficiencies,  the  EMS  is  responsible  for  finding  the  most 
appropriate  power  splitting  of  Pu.  This  power  splitting  process  will 
lead  to  the  generation  of  the  power  control  signals  Pfc  (for  the  fuel 
cell  power  system)  and  P£,  (for  the  Propulsion  System).  Given  these 
power  control  signals,  the  uni-directional  DC— DC  converter  will 
provide  the  corresponding  electric  power  Pfc  and  the  Propulsion 
System  will  deliver  the  corresponding  Pm.  Pb  represents  the  battery 
power. 

This  architecture  allows  the  battery  power  source  to  be  fully 
responsible  to  track  the  driver  power  demand.  At  the  same  time, 
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the  fuel  cell  operates  at  its  maximum  efficient  in  order  to  save 
hydrogen.  Thus,  the  fuel  cell  power  production  management 
problem  is  to  find  the  start  and  stop  time  sequence  during  the 
vehicle  motion  in  order  to  minimize  the  hydrogen  consumption 
and  to  keep  the  battery  pack  in  a  healthy  state.  In  the  paper,  we 
assume  that  the  battery  healthy  state  is  negatively  affected  by  the 
depth-of-discharge  (DOD)  and  by  the  high-power  discharge  rates: 
high  DOD  and  high  discharge  rates  will  shorten  the  battery  life 
span.  Clearly,  the  role  of  the  PEMFC  in  this  architecture  is  to  support 
the  battery  pack  by: 

•  reducing  the  DOD; 

•  reducing  the  battery  discharge  rate; 

•  preventing  the  battery  to  be  depleted  below  a  minimum 
prescribed  energy. 


flow  rate,  respectively.  The  fuel  cell  power  system  efficiency  rjfc  is 
defined  as: 


Vfc(k)  = 


Pfc(fe) 

Kmrhm(k) 


(7) 


One  of  the  goals  of  the  Fuel  Cell  is  to  reduce  the  impact  of  high 
Battery  Pack  discharge  rate  due  to  the  propulsion  system  fast 
dynamics.  These  discharge  periods  are  mostly  related  to  the  fast 
dynamics  of  the  user  power  demand  Pu,  compared  to  the  Fuel  Cell 
dynamics.  So  the  fuel  cell  power  source  dynamics  need  to  be 
included  into  the  Energy  Management  System.  Equation  (8) 
represents  a  first-order  model  of  a  PEMFC  [20]: 


Pfc(k  +  1)  =A(k)Pfc(k)  +  B(k)P*fc(k) 


(8) 


3.  FC-PHEV  longitudinal  model  for  energy  management 

3.1.  Energetic  model 


where  0  <  A(k)  <  %  B(k)  and  Pfc(k)  represent  respectively  the  Fuel 
Cell  dynamic  constant,  the  input  coefficient  and  the  input  power 
command  (see  Fig.  1).  The  sequences  of  A  and  B  can  be  determined 
through  any  classical  system  identification  methods. 


The  longitudinal  motion  of  the  Fuel  Cell  PHEV  is  represented  by 
Equation  (1)  [22]  and  the  corresponding  mechanical  power  is 
represented  by  (2). 

Fv(k)  =  Mvv(k)  +~paCdAvv^(k)  +  Mvgsin(S(k))  +  Mvgficos(d(k)) 

(1) 

where  k,  Mu,  g,  pa,  p,  Cd  and  Au  represent  respectively  the  sampling 
index,  the  vehicle’s  total  mass,  the  gravity  constant,  the  air  density, 
the  rolling  resistance  coefficient,  the  air  drag  coefficient  and  the 
active  frontal  area;  6,  i)  and  Ud  represent  respectively  the  road  grade, 
the  longitudinal  acceleration  and  the  vehicle  speed  relative  to  the 
air.  0  <  k  <  N  where  N  is  the  number  of  samples  in  a  given  driving 

Pm(k)  =  Fv(k)v(k)  (2) 

The  electric  power  Pe  required  for  the  vehicle  propulsion 
depends  mainly  on  the  Electric  Propulsion  System  efficiency  ?jm  and 
Pm- 

Pe(k)  =  VmHk)Pm(k)  (3) 


3.2.  Fuel  Cell  and  battery  pack  efficiency  models 

The  Fuel  Cell  and  the  Battery  Pack  efficiency  degradations  are 
currently  being  investigated  in  the  research  community  [27—30], 
However,  the  reported  results  demonstrate  that  these  phenomena 
are  difficult  to  model  precisely.  So,  to  analyze  the  impact  of  the 
efficiencies  on  the  overall  management  system,  we  use  a  semi- 
empirical  approach  based  on  the  available  power  sources  manu¬ 
facture  data. 

3.2.1.  Fuel  Cell  efficiency 

The  Fuel  Cell  degradation  depends  on  several  factors: 
number  of  voltage  cycling,  power  level,  power  dynamics  and 
ambient  temperature  and  operating  temperature.  However,  the 
detail  interaction  of  these  factors  with  the  degration  is  complex. 
So  we  assume  that  for  a  constant  power  production,  the  Fuel 
Cell  efficiency  will  decayed  exponentially  with  the  time  [21,31] 
as  represented  by  Equation  (9).  Even  if  the  variable  level  of 
the  FC  power  is  used,  a  detained  relationship  regarding  the 
power  level,  the  efficiency  and  the  degradation  can  be  used  as 
well. 


The  Equation  (4)  represents  the  power  balance  at  the  electric 
node  (see  Fig.  1).  A  positive  value  of  the  bi-directional  converter 
power  Pb  represents  the  power  from  the  Battery  Pack  whereas 
a  negative  value  of  Pb  represents  the  power  to  the  Battery  Pack. 

Pb(k)  =  rm'(k)Pm(k)-Pfc(k)  (4) 

The  relationship  between  Eb  and  Pb  is  represented  the  following 
equation: 

Eb(k+  1)  =  Eb(k)  -  r,b\k)Pb(k)\T  (5) 

where  k,  AT  and  Tjb  represent  the  sampling  index,  the  sample  time 
period  and  the  battery  power  source  efficiency  (this  efficiency 
includes  the  Bi-directional  Converter  and  the  Battery  Pack  effi¬ 
ciencies),  respectively. 

On  an  energetic  point  of  view  and  assuming  that  the  fuel  cell 
parameters  are  adequately  controlled,  the  following  equation 
represents  the  energy  balance  in  the  Hydrogen  Tank  [20,23,24]: 

Efc(k  +  1)  =  Efc(k)  -  Kmmm(k)\T  (6) 

where  £fC,  Km  and  mm  represent  the  hydrogen  energy  in  the 
storage,  the  hydrogen  high  heating  value  [25,26]  and  the  hydrogen 


Vtc (k)  =  €e  D^T  (9) 

where  rft.  and  Dfc  represent  the  fuel  cell  power  source  initial  effi¬ 
ciency  and  the  efficiency  decay  rate,  respectively.  AT  is  the  sampling 
time. 

We  consider  that  the  manufacturer  provides  the  Fuel  Cell  life¬ 
time  tfc.  Knowing  the  Fuel  cell  initial  efficiency  (i.e.  the  measured 
efficiency  when  the  Fuel  Cell  was  used  for  the  first  time),  Dfc  is 
obtained  using  the  following  equation: 


where  rjfc  represents  the  minimum  value  corresponding  to  the  Fuel 
Cell  end-of-life  efficiency. 


3.2.2.  Battery  pack  efficiency 

The  battery  pack  efficiency  depends  on  several  factors.  The 
discharge  power  (or  the  discharge  current)  and  the  depth  of 
discharge  (DOD)  are  among  the  most  influential  parameters  of  the 
Battery  Pack  efficiency  degradation  [32,33].  Often,  the  manufacture 
provides  two  different  maps  [34]: 
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updated  at  the  end  of  the  driving  cycle  i  with  the  following 
equation  and  the  battery  is  recharged  completely  for  the  next 
driving  cycle: 


DOD(i)  = 


£"=0fb(fc)AT 

Cb(0) 


N  is  the  number  of  current  samples  for  a  given  driving  cycle  i. 

For  each  AT,  the  instantaneous  efficiency  ?j£(k)  is  retrieved  using 
the  map  Fj.  However,  the  value  of  ?jb(fc)  does  not  take  into  account 
the  battery  degradation  effect  represented  by  the  map  F2.  To 
include  this  effect,  we  define  the  degradation  coefficient  as  follows: 


where  tb  (t  —  1)  is  the  remaining  lifetime  after  the  previous  driving 
cycle.  0  <  r)£(k)  <  I,  vk. 

The  Battery  Pack  efficiency  is  obtained  with  the  following 
equation: 


Fig.  3.  Experimental  Fuel  Cell  startup  delay  during  a  high  discharge  phase:  Pu  and  Pfc 
represent  the  power  demand  and  the  fuel  cell  power,  respectively. 

•  the  battery  capacity  Cb  as  a  function  of  the  discharge  current 
Ib(k).  Knowing  the  battery  fully  charge  capacity  Cb(0),  we 
defined  the  energetic  efficiency  map  rjg  =  F1(Cb/Cb(0))  ; 

•  the  expected  lifetime  (or  expected  number  of  cycles)  as 
a  function  of  the  DOD:  tg  =  F2(DOD(i));  i  represents  the 
battery  charge/discharge  cycle  index.  The  value  of  DOD(i)  is 


rib(k)  =  77g(k)r7g(k)  (10) 

Substituting  P^k)  and  ? jb(k)  in  Equation  (5)  with  Equations  (4) 
and  (10),  respectively,  yields: 

Eb(k+l)=£b(k)—  [t/b(k)t;b(k)]  1  [Vm\k)Pm(k)-Pfc(k)\  AT  (11) 

Taking  into  account  the  Fuel  Cell  efficiency,  the  hydrogen  energy 
in  the  tank  is  always  given  by  the  following  equation: 


Time(s) 

Fig.  4.  Three  different  synthetic  driving  cycles:  (a)  US06  Driving  Cycle,  (b)  New  European  Driving  Cycle  and  (c)  Japanese  10/15-Mode  Driving  Cycle. 
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Fig.  5.  Three  acceleration  profiles:  (a)  US06  Driving  Cycle,  (b)  New  European  Driving  Cycle  and  (c)  Japanese  10/15-Mode  Driving  Cycle. 


E{c(k  + 1)  =  EfC(k)  -  [4e-D^r]  -1PfcAT  (12) 

Clearly,  the  hydrogen  consumption  and  the  Battery  Pack  energy 
consumption  increase  as  the  different  power  source  efficiencies 
decrease  over  a  long  period  of  time.  Since  the  amount  of  energy  is 
limited,  the  following  constraints  must  hold: 

EfT  <  Eft(k)  <  4“  (13) 


Discharge  Current  (A) 

Fig.  7.  A  lead-acid  battery  efficiency  curve  [46], 
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Fig.  8.  Battery  lifetime  map  [46]. 


Cn  <  Eb(k)  <  erx  (W) 

where  £™n,£™ax,£“m  and  Eg1™  represent  the  minimum  hydrogen 
energy  in  the  tank,  the  maximum  hydrogen  stored,  the  minimum 
Battery  Pack  energy  and  the  maximum  Battery  Pack  energy.  Eg™  is 
defined  as  the  minimum  energy  threshold  that  should  be  kept  in 
the  battery  Pack  in  order  to  reduce  its  degradation.  This  value  will 
play  an  important  role  for  the  Energy  Management  System. 

4.  Energy  management  problem  formulation 

4.1.  High  power  discharge  anticipation  problem 

One  of  the  goals  of  the  fuel  cell  power  source  is  to  reduce  the 
impact  of  high  discharge  on  the  Battery  Pack.  The  high  discharge 
usually  occurs  during  the  vehicle  acceleration  phase.  If  the  Fuel  cell 
dynamics  are  similar  (i.e.  as  quick  as)  to  the  vehicle  dynamics,  then 
it  would  be  easier  to  start  it  at  the  beginning  of  these  acceleration 
phases.  Unfortunately,  the  slow  Fuel  Cell  startup  time  imposes  to 
anticipate  the  upcoming  acceleration  phases.  Fig.  3  illustrates 
a  typical  high  discharge  during  the  vehicle  acceleration.  The  fuel 
cell  power  reaches  its  maximum  several  second  after  the  end  of  the 
high  discharge  phase.  This  example  shows  that  the  Fuel  Cell  should 
be  started  several  seconds  before  this  high  discharge  phase  begins 
in  order  to  reduce  the  negative  impact  on  the  Battery  Pack. 

4.2.  Real-time  power  splitting  problem 

In  motor  mode,  the  Electric  Propulsion  System  is  supplied  with 
the  electric  power  Pe  from  the  on-board  power  sources  (Fuel  Cell  or 
Battery  Pack:  the  arrows  with  the  solid  lines  in  Fig.  1).  At  each 
timestamp  k,  the  human  driver  provides  the  desired  power 
command  Pu(k).  Given  Pu(k),  the  Energy  Management  System  will 
generate  the  required  optimal  control  signals  for  the  fuel  cell  power 
source  Pfc(k)  and  for  the  Propulsion  System  P^fk).  We  propose  in 

Table  1 

PHEV  parameters. 

Parameter  Value  Parameter  Value 

~  1200  kg  A,  1.2  m2 

pa  1.2  kgm  3  p  0.018 

Cd  0.5 


Table  2 

Simulation  parameters. 


Parameter 

Value 

Parameter 

Value 

Eb(0) 

100  Wh 

Ru 

100/36002 

Vm 

0.9 

Rfc 

10/36002 

£mi" 

50  Wh 

C 

1 

£max 

100  Wh 

4 

20  Ws 

;Lr. 

100  Wh 

Vfc 

0.5 

1*111 

0 

Dfc 

10~9 

prnax 

1500  W 

AT 

1  s 

A[k) 

0.72,  Vfc 

B{k) 

0.95,  Vk 

this  section,  a  problem  formulation  that  will  allow  a  real-time 
generation  of  these  optimal  control  signals  in  order  to  fulfill  at 
least  to  the  following  conditions. 

•  The  Battery  Pack  energy  Eb  should  be  kept  above  the  safe 
energy  threshold  £™n.  Complying  with  this  condition  will  help 
to  reduce  Battery  Pack  degradation  [35,36], 

•  The  fuel  cell  power  source  should  contribute  to  reduce  the  high 
power  discharges  from  the  Battery  Pack.  Complying  with  this 
condition  will  have  two  benefits:  the  improvement  to  the 
battery  power  source  efficiency  and  the  reduction  of  the 
Battery  Pack  degradation  [37], 

•  The  Fuel  Cell  must  operate  at  its  maximum  efficiency  in  order 
to  save  hydrogen. 

One  way  to  simultaneously  take  into  account  the  above 
mentioned  conditions  is  to  formulate  an  optimization  problem 
with  a  multicriteria  cost  function.  The  proposed  cost  function  (see 
Equation  (15))  has  two  components.  The  first  component  L(k ) 
(Equation  (16))  is  the  immediate  cost  whereas  the  second  compo¬ 
nent  i p  (Equation  (17))  represents  the  cost  that  will  result  in  the 
application  of  the  selected  control  signals  P^(k)  and  Pfc(k). 

J(k)  =  L(k)  +  <f>  (15) 

where: 


Lik)  =  1  ((Eb(k)  -  E^y+Ru(Pu(k)  -  Pm(k))2+Rfc(pfc(k) 


Fig.  9.  Experimental  small  utility  truck  FC-HEV. 


412 


S.  Kelouwani  et  al.  /  Journal  of  Power  Sources  221  (2013)  406-418 


and  where: 

0  =  i(£b(k+l)-£™n)2  (17) 

and  where  Ru  >  0  is  the  penalty  coefficient  of  the  deviation  between 
the  driver  power  command  and  the  actual  power  command  sent  to 
the  propulsion  system;  RfC  >  0  is  the  penalty  coefficient  of  the 
deviation  between  the  fuel  cell  power  and  the  anticipated  power  Pa 
based  on  the  upcoming  high  discharge  from  the  Battery  Pack;  C  is 
the  result  of  the  division  between  the  energy  costs  of  hydrogen  and 
electricity.  We  assume  that  the  hydrogen  energy  is  more  expensive 
than  the  electricity  energy.  C  >  1  is  the  hydrogen/electricity  cost 
ratio  applied  to  the  hydrogen  consumption  mH2 •  According  to 
Equations  (7)  and  (9),  the  expression  of  mH2  is: 

ifefsS*  (,8) 

From  Equation  (16),  the  minimum  of  L(k)  is  obtained  when: 

•  the  Battery  Pack  energy  is  close  to  the  safe  minimum  value. 
This  criterion  will  allow  the  Energy  Management  System  to 
operate  in  charge  depleting  mode  until  £b  reaches  £[™n:  in 
addition,  since  C  >  1,  the  battery  energy  will  be  preferably  used 
instead  of  the  hydrogen  energy 

•  the  driver  power  demand  Pu  is  close  to  the  mechanical  power 
Pm  provided  to  the  vehicle  dynamics;  this  criterion  will  permit 
the  vehicle  to  behave  dynamically  as  desired  by  the  driver; 

•  the  hydrogen  consumption  rate  is  very  low;  this  criterion  will 
allow  the  use  of  the  minimum  hydrogen  in  order  to  fulfill  to  the 
other  criteria; 

•  the  fuel  cell  power  Pfc  is  close  to  the  anticipated  high  demand 
from  the  driver  Pa;  as  said  before,  the  dynamics  of  the  Fuel  Cell 


is  slower  than  the  power  demand  dynamics.  So,  it  is  necessary 
to  anticipate  the  next  upcoming  high  demand  in  order  to  allow 
enough  time  to  the  Fuel  Cell  to  be  powered  up  and  contribute 
to  the  reduction  of  this  high  power  discharge.  A  method  for 
computing  Pa  is  given  in  the  next  section. 

The  second  component  of  J(k)  represented  by  <p  (Equation  (17)) 
is  added  in  order  to  avoid  depleting  to  much  the  Battery  Pack,  when 
the  £b  becomes  close  to  £[)lln. 

Hence  the  optimization  problem  is  formulated  as  follows:  for 
each  timestamp  k,  find  the  control  signals  Pfc(k)  and  P^fk)  so  that 
the  cost  function  (16)  is  minimized  under  the  constraints  repre¬ 
sented  by  the  following  expressions: 

1.  £b(k+  1)  =  £b(k)  -  [r;ge-D^T]-1[77-1(k)P^(k)  -  Pfc(k)]AT 

2.  £fc(k  + 1)  =  EfC(k)  -  [7/fCe_D(cfcAT] -1  PfC AT 

3.  Pfc(k  + 1)  =  A(k)Pfc(k)  +  £(k)P;c(k) 

4.  E™in  <  Efc(k)  <  £™x 

5.  £™n  <  £b(k)  <  £™ax 

6.  In  order  to  let  the  fuel  cell  to  operate  at  its  maximum  efficiency, 
the  value  of  Pfc  must  belong  to  (0,  Pj™x }  where  P™x  represents 
the  power  produced  with  the  maximum  efficiency. 

7.  It  is  desirable  to  limit  the  number  of  the  Fuel  Cell  start/stop 
cycles  in  order  to  reduce  its  degradation  [38].  So  the  energy 
produced  by  the  Fuel  Cell  during  each  start/stop  cycle  should  be 
greater  than  a  minimum  value  E*fc  [30].  A  large  value  of  E*fc  will 
let  the  Fuel  Cell  to  work  during  a  long  period.  As  a  result,  the 
number  of  start/stop  cycles  will  be  smaller  than  if  the  E/c  is 
small. 

In  the  next  section,  we  propose  a  new  method  to  solve  this 
optimization  problem. 
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5.  Optimal  power  splitting  solution 

5.1.  High  power  discharge  anticipation  algorithm 

The  first  step  for  solving  the  aforementioned  optimization 
problem  is  to  determine  the  sequence  of  anticipated  high  power 
discharge  Pa(k),  k  —  1,2,3,...  One  way  to  find  Pa  is  to  train  a  neural 
network  in  order  to  predict  the  next  occurrence  of  high  power 
discharge.  However,  the  convergence  rate  may  be  very  slow.  Thus, 
we  propose  a  simple  and  efficient  method  to  find  Pa. 

The  high  power  discharge  from  Battery  Pack  is  mainly  due  to  the 
vehicle  acceleration  imposed  by  the  driver.  Often,  this  acceleration 
happens  in  two  situations.  The  first  one  is  related  to  the  vehicle 
motion  transition  from  a  deceleration  phase  to  an  acceleration 
phase.  In  this  case,  the  vehicle  is  doing  a  stop-and-go  or  a  slow- 
down-and-go.  The  second  situation  occurs  when  the  vehicle 
accelerates  from  a  steady  speed.  In  practice,  the  first  situation 
occurs  most  of  the  time.  So,  it  is  possible  to  anticipate  the  upcoming 
acceleration  by  observing  when  the  deceleration  occurred. 

Fig.  4  shows  three  commonly  used  synthetic  driving  cycles:  the 
high  acceleration  aggressive  driving  schedule  of  the  United  States 
Environmental  Protection  Agency  (US06)  [39,40],  the  New  Euro¬ 
pean  Driving  Cycle  (NEDC)  [41]  and  the  Japanese  10/15  Mode 
Driving  Cycle  Q1015)  [42],  The  acceleration  profiles  obtained  from 
these  driving  cycles  are  shown  in  Fig.  5  and  suggest  that  most  of  the 
time,  a  deceleration  phase  is  followed  by  an  acceleration  phase 
until  the  cycles  are  completed. 

Hence,  Fig.  6  represents  the  algorithm  used  to  generate  the 
anticipate  high  power  discharge  signal  Pa  based  upon  the 


knowledge  about  the  driver  power  demand  Pu.  We  consider  that 
a  deceleration  phase  happens  when  Pu  <  0.  It  is  important  to 
mention  that  the  algorithm  is  not  trying  to  predict  when  the  driver 
will  stop.  Instead,  it  observes  the  deceleration  phase  and  uses  this 
knowledge  to  trigger  the  fuel  cell  operation.  One  can  start  the  Fuel 
Cell  at  the  beginning  of  a  deceleration  phase  (when  Pu  <  0). 
However,  the  positive  value  of  the  Fuel  Cell  power  could  limit  the 
ability  of  the  battery  to  properly  handle  the  regenerative  power 
that  may  be  produced  during  the  vehicle  deceleration.  To  avoid  this 
drawback,  we  propose  to  start  the  Fuel  Cell  at  the  end  of  the 
deceleration  phase  which  is  represented  by  the  transition  between 
a  negative  value  of  Pu  and  the  following  positive  value  of  Pu.  In 
operation,  the  stack  must  be  set  to  produce  the  power  at  its 
maximum  efficiency  P™ax.  The  Fuel  Cell  is  stopped  at  the  end  of  the 
acceleration  phase. 

The  anticipatory  algorithm  assumes  that  if  the  difference 
between  the  current  high-power  demand  value  and  the  previous 
one  is  less  than  a  given  threshold  rp  <  0  when  Pu  >  0,  then  the  high 
demand  is  decreasing.  So,  the  Fuel  Cell  can  be  stopped.  One  possible 
way  to  determine  rp  <  0  is  to  consider  the  average  power  demand 
and  use  the  following  equation: 


=  Ej=opuO)  , 


(19) 


The  proposed  algorithm  is  real-time  because  no  a  prior  knowl¬ 
edge  of  the  entire  driving  cycle  is  required  and  all  involve  pro¬ 
cessing  are  simple  algebraic  operations.  Given  the  sequence  of 

Pu(k),  k  —  1,  2 . N,  the  optimal  solution  of  the  splitting  problem  is 

presented  in  the  next  section. 
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5.2.  Locally  optimal  power  splitting  solution 

Several  methods  exist  to  solve  this  static  and  nonlinear  opti¬ 
mization  problem  [43]:  iterative  descent  methods,  approximation 
methods  and  dual  and  primal-dual  methods.  We  have  used  the  line 
search  algorithm  for  nonlinear  programming  as  outlined  in 
Ref.  [43],  This  algorithm  is  selected  because  it  provides  an  iterative 
optimization  method  that  is  simple  and  efficient  to  handle  for  real¬ 
time  applications. 

Line  search  method:  The  generic  algorithm  of  line  search  is  an 
iterative  method  which  involves  two  main  steps  [44]: 

•  the  search  of  a  descent  direction  along  the  cost  function  J(k); 

•  the  selection  of  a  step  size  along  the  selected  direction. 

In  this  paper,  we  used  an  enhanced  line  search  algorithm  orig¬ 
inally  proposed  in  Ref.  [45].  T  delimits  the  control  signal  search 
domain  on  which  the  optimization  is  performed.  This  domain  is 
divided  is  (Nfc  x  Nm)  control  signal  grid  points  where  Nfc  and  Nm 
represent  the  number  of  equal  intervals  on  the  PfC  dimension  and 
on  the  Pm  dimension,  respectively.  Using  each  grid  point  as  an 
initial  guess,  the  line  search  algorithm  described  in  Ref.  [45]  is  used 
in  order  to  find  local  minimums.  The  optimal  control  signals  P'fc  and 
P*m  are  control  signals  Pfc  and  Pm  that  produce  the  smallest  value  of 
J(k)  by  considering  all  minimums  found  using  the  grid  points  as  the 
initial  guesses  and  by  considering  all  constraints  involved  in  the 
problem  formulation. 

The  controller  obtained  uses  2000  grid  points  (Nfc  =  2, 
Nm  —  1000).  This  number  of  points  allows  a  real-time  search  of  the 
optimums  Pfc  and  Pj;,.  For  all  performed  simulation,  the  maximum 
execution  time  of  the  line  search  method  is  0.98  ms  on  a  Personal 
Computer  (PC)  with  Pentium  IV  3  GHz. 


6.  Simulation  and  discussion 

The  simulation  pursues  three  main  objectives:  (i)  the  validation 
of  the  Fuel  Cell  power  production  through  the  anticipative  algo¬ 
rithm;  (ii)  the  analysis  of  the  impact  of  the  power  system  efficiency 
deteriorations  on  the  overall  real-time  Energy  Management  System 
proposed  in  this  paper;  (iii)  the  comparative  study  of  the  energy 
saving  between  the  proposed  method  and  a  typical  rule-based 
energy  management  method  [2], 

6.1.  Simulation  context 

The  vehicle  and  the  simulation  parameters  are  summarized  in 
Table  1  and  Table  2,  respectively.  These  vehicle  parameters  corre¬ 
spond  to  the  small  utility  truck  shown  in  Fig.  9  and  which  is  being 
used  in  our  research.  We  used  a  scale  down  factor  of  0.01  on  the 
power  request  and  the  different  energy  storage  in  order  to  match  the 
fuel  cell  rated  power  of  Fig.  9.  The  vehicle  has  an  on-board  Fuel  Cell 
system  (with  a  rated  power  of  1500  W)  to  recharge  the  Battery  Pack. 
A  72  V  lead-acid  battery  pack  is  used  as  the  electric  energy  storage. 
For  the  simulation,  we  considered  that  the  FC  response  delay  is  5  s.  At 
the  beginning  of  each  simulation,  we  assume  that  the  Battery  Pack  is 
fully  charged  and  no  a  prior  knowledge  of  the  driving  cycle  is 
available.  The  efficiency  map  and  the  lifetime  map  are  represented  in 
Figs.  7  and  8,  respectively  [46,47],  The  system  must  rely  only  on  the 
current  and  the  past  data  collected.  The  Fuel  Cell  system  is  also 
assumed  to  be  properly  controlled  by  a  PID  controller. 

6.2.  Fuel  Cell  anticipative  power  production  validation 

To  validate  the  Fuel  Cell  power  production  with  the  proposed 
anticipative  algorithm,  we  consider  the  driving  cycle  represented  in 
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Fig.  10  (graph  (a))  and  which  is  obtained  with  the  vehicle  shown  in 
Fig.  9.  This  driving  cycle  has  six  stops-and-go,  which  reflect  a  typical 
usage  of  such  a  utility  truck.  The  blue  curve  (in  web  version)  in  the 
graph  (b)  (Fig.  10)  represents  the  power  demand  Pu.  Except  for  the 
starting  high-power  demand  (timestamp  k  —  0),  where  the  Fuel  Cell 
was  not  powered  on,  we  observe  in  the  graph  (b)  of  Fig.  10  that  the 
anticipative  algorithm  can  track  each  subsequent  high-power 
demand  from  the  driver.  Furthermore,  the  Fuel  cell  starts  at  the 
end  of  each  deceleration  phase  and  stops  as  soon  as  the  corre¬ 
sponding  high-power  demand  is  ending.  As  a  result,  the  peak 
power  drawn  from  the  Battery  Pack  is  reduced  by  29.5%  in  average 
(see  the  green  curve  (in  web  version)  in  the  graph  (b)).  In  addition, 
the  Fuel  Cell  power  production  does  not  interfere  with  the  regen¬ 
erative  power  (the  negative  part  of  Pb  coincides  with  the  negative 
part  of  Pu).  Due  to  its  slow  dynamics,  if  two  stops-and-go  are  too 
close  to  each  other,  the  Fuel  Cell  will  be  powered  on  at  the  begin¬ 
ning  of  the  first  one  and  will  be  stopped  at  the  end  of  the  decel¬ 
eration  phase  of  the  second  one.  This  fuel  cell  control  behavior  is 
observed  during  the  time  interval  between  155  s  and  168  s  (see  the 
graph  (b),  Fig.  10). 

6.3.  Impact  of  the  hydrogen/electricity  cost  ratio  and  the  power 
source  efficiency  deterioration 

Two  important  parameters  influence  the  energy  management 
system  and  the  power  splitting  algorithm:  the  power  source  effi¬ 
ciencies  and  the  hydrogen/electricity  cost  ratio.  To  analyze  the 
impact  of  the  battery  power  system  efficiency  rjb,  we  used  the 


driving  cycle  presented  in  Fig.  10  (in  the  graph  (a)).  We  conducted 
two  different  simulations  in  order  to  demonstrate  the  effectiveness 
of  the  proposed  energy  management  and  power  splitting  methods: 

•  in  the  first  simulation,  we  assume  that  the  Battery  Pack  is  new. 
Thus,  the  efficiency  r/b  is  high; 

•  in  the  second  simulation,  the  Battery  Pack  degradation  is  high 
and  the  efficiency  is  low. 

During  simulation,  Tjfc  is  assumed  to  change  slowly  (see 
parameters  in  Table  2). 

6.3.1.  Battery  power  system  with  high  efficiency 

Fig.  11  shows  the  Fuel  Cell  power  generation  when  the  two 
efficiencies  are  high.  We  see  in  the  graph  (c)  of  this  figure  that  the 
Battery  Pack  energy  is  most  of  the  time  greater  or  close  to  the 
minimum  energy  to  keep  in  the  batteries  (Eg1'"  =  50  Wh).  In 
addition,  the  Fuel  Cell  can  compensate  for  the  high-power  demand 
due  to  the  stop-and-go  as  shown  in  the  graph  (b)  of  Fig.  11. 

The  cost  function  (15)  has  two  conflicting  requirements:  (i) 
maintain  the  battery  energy  greater  than  Eg™  by  using  the  fuel  cell 
energy;  (ii)  minimize  the  hydrogen  energy  cost  which  is  repre¬ 
sented  (C  x  jt1h2)-  Note  that  C  >  1  is  set  by  the  user  as  an  input 
parameter.  When  the  hydrogen/electricity  cost  ratio  C  increases, 
the  Energy  Management  System  reduces  the  use  of  hydrogen  by 
decrease  the  Fuel  Cell  power  production  in  order  to  minimize  the 
cost  function  represented  by  Equation  (15).  Flence,  the  Battery  Pack 
is  depleted  in  order  to  meet  the  power  demand.  There  exist  a  cost 
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Fig.  13.  Comparative  study  with  US06  Driving  Cycle:  (a)  power  demand  Pu,  Fuel  Cell  optimal  power  production  Pfc  and  Fuel  Cell  power  production  with  the  rule-based  method  P^f ; 
(b)  battery  pack  optimal  energy  profile  Ejj  and  battery  pack  rule-based  energy  profile  Eg6 ;  (c)  stored  hydrogen  with  the  optimal  energy  profile  Ejc  and  stored  hydrogen  with  the  rule- 
based  energy  profile 
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value  for  which  the  Fuel  Cell  contribution  cannot  help  to  maintain 
the  Battery  Pack  final  energy  £b(JV)  close  to  £™n.  We  designated  this 
specific  cost  value  as  the  critical  cost.  By  increasing  the  cost  ratio 
from  1  to  100,  we  found  that  this  critical  cost  ratio  is  10  for  the 
power  system  with  high  efficiencies.  Thus,  if  the  cost  ratio  is  lower 
than  10,  the  Energy  Management  System  can  maintain  the  healthy 
level  of  the  Battery  Pack  using  the  Fuel  Cell  power  production. 

6.3.2.  Battery  power  system  with  low  efficiency 

When  the  battery  is  aging,  it  efficiency  decreases.  So  more 
hydrogen  is  needed  to  maintain  its  minimum  energy  above  Eg"n 
during  the  driving  cycle  and  the  Energy  Management  System  must 
adapt  its  power  splitting  algorithm  accordingly.  Fig.  12  shows 
a  typical  result  with  a  low  efficiency  for  the  battery.  The  operating 
time  between  a  start  and  a  stop  of  the  Fuel  Cell  is  longer  when  the 
Battery  Pack  efficiency  is  low  compared  to  the  previous  case  for 
which  the  efficiency  is  high  (see  the  graph  (b)  in  Fig.  12  and  the 
graph  (b)  in  Fig.  11 ).  In  addition,  the  battery  energy  is  kept  greater 
or  close  to  Eg™  during  the  driving  cycle. 

To  analyze  the  effect  of  the  hydrogen/electricity  cost  ratio  on  the 
capacity  to  maintain  the  battery  depletion  close  to  £™n,  we  run 
several  simulations  with  different  values  of  C.  These  results  indicate 
that  the  value  of  the  stored  hydrogen  at  the  end  of  the  driving  cycle 
Efc(N)  are  lower  than  their  corresponding  values  obtained  when  the 
battery  efficiency  is  high.  Indeed,  as  the  battery  power  source  is  less 
efficient  than  the  Fuel  Cell,  more  hydrogen  is  used  to  produce  the 
propulsion  energy.  Thus,  the  Energy  Management  System  can 
reduce  the  use  of  the  less  efficient  power  source  in  favor  of  the  most 
efficient  one.  Since  the  Battery  Pack  energy  is  reduced  in  the  low 


efficiency  case,  the  cost  ratio  corresponding  to  the  critical  point 
(C  =  36)  is  much  greater  than  the  cost  ratio  corresponding  to  the 
critical  point  obtained  when  the  efficiency  is  high. 

6.4.  Comparative  study 

We  consider  the  high  efficiency  for  the  Battery  Pack  and  three 
commonly  used  driving  cycles  (see  Fig.  4).  The  aim  of  the 
comparative  study  is  to  show  the  effectiveness  of  the  proposed 
Energy  Management  System  by  analyzing  the  Eb  and  Efc  dynamics 
when  one  of  the  three  driving  cycles  is  used.  For  comparative 
purposes,  a  commonly  used  rule-based  method  for  real-time 
energy  management  is  considered.  This  method  lets  the  Fuel  Cell 
to  start  and  to  stop  if  Eb  is  lower  than  a  predefined  lower  limit  and  if 
Eb  is  higher  than  a  predefined  high  limit,  respectively.  For  the 
comparative  tests,  the  upper  limit  of  Eb  is  set  to  Ejjp  =  55  Wh  and 
lower  limit  is  set  to  e£ov™  =  45  Wh.  EjJ1™  is  set  to  50  Wh.  Since  the 
Fuel  Cell  maximum  power  is  limited  to  1500  W  and  the  average 
power  demand  using  the  different  driving  cycles  is  much  greater, 
we  introduced  a  scale  down  factor  of  0.01. 

Figs.  13-15  show  the  comparative  study  results.  The  hydrogen 
energy  used  for  the  three  driving  cycles  (graph  (c)  of  Figs.  13—15) 
when  the  optimization  method  proposed  in  this  paper  used  is 
less  than  the  ones  obtained  with  the  rule-based  method. 
Furthermore,  for  the  three  driving  cycles,  the  Eb  at  the  end  of  the 
cycles  is  always  greater  or  close  to  £™n  when  the  optimization 
method  is  used.  These  results  demonstrate  that  the  proposed 
method  is  efficient  in  preventing  the  Battery  Pack  to  be  over 
depleted. 
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Fig.  15.  Comparative  study  with  J10/15  Driving  Cycle:  (a)  power  demand  Pu,  Fuel  Cell  optimal  power  production  Pj.  and  Fuel  Cell  power  production  with  the  rule-based  method  P{® ; 
(b)  battery  pack  optimal  energy  profile  E^  and  battery  pack  rule-based  energy  profile  Ef6 ;  (c)  stored  hydrogen  with  the  optimal  energy  profile  Ej.  and  stored  hydrogen  with  the  rule- 
based  energy  profile  Eji? 


Degradation  analysis  i 


Efficiency  US06  NEDC  J10/15 


Optimized  Rule-based  Optimized  Rule-based  Optimized  Rule-based 

rrf  after  100  h  98.6%  98.0%  99.3%  99.0%  97.0%  98.3% 

Tjjj  after  200  h  96.7%  95.7%  98.4%  97.4%  96.5%  94.4% 


To  analyze  the  effect  of  the  Battery  Pack  degradation  on  the 
proposed  energy  management  system,  we  ran  three  simulations 
using  the  extended  versions  of  the  three  driving  cycles  (US06, 
NEDC,  J10/15).  Each  extended  driving  cycle  is  obtained  by  concat¬ 
enating  the  original  driving  cycle  in  order  to  obtain  a  total  of 
5000  h.  Using  the  battery  efficiency  curve  (Fig.  7)  and  the  battery 
lifetime  map  (Fig.  8),  we  investigate  the  degradation  when  the  rule- 
based  algorithm  or  the  optimal  method  is  used.  We  considered  that 
when  the  DOD  is  greater  than  75%,  the  Battery  Pack  is  completely 
recharged  and  the  hydrogen  tank  if  fully  filled.  The  following  table 
summarizes  the  results. 

In  Table  3,  values  when  the  optimized  method  is  used  is 
greater  than  the  corresponding  values  when  the  rule-based 
method  is  used.  These  results  suggest  that  the  optimal  energy 
management  system  is  more  effective  to  preserve  the  Battery  Pack 
lifetime  than  the  rule-based  method. 

7.  Conclusion 

A  strategy  for  a  Fuel  Cell-PHEV  based  on  optimal  power  splitting 
between  the  Battery  Pack  and  the  Fuel  Cell  is  proposed  and 


discussed.  This  power  splitting  method  aims  at:  (i)  reducing  the 
high-power  demand  on  the  batteries,  especially  during  the  vehicle 
stop-and-go  operations  and  (ii)  preventing  the  excessive  depletion 
which  may  shorter  the  Battery  Pack  life  span.  An  experimental 
small  utility  truck  longitudinal  model  was  built  and  used  for  the 
control  design  and  analysis.  To  reduce  high-power  demand,  the 
Energy  Management  System  uses  an  anticipatory  control  tech¬ 
nique.  Furthermore,  a  multi-criteria  optimization  including  the 
hydrogen/electricity  cost  ratio  and  the  anticipatory  control  is 
formulated  and  solve  using  the  line  search  method.  The  overall 
control  system  does  not  require  a  prior  knowledge  of  the  trip 
information.  Therefore,  the  real-time  implementation  of  this 
Energy  Management  System  is  realistic.  The  analysis  of  the  effects 
of  the  efficiency  degradations  on  the  system  demonstrated  that 
there  exist  a  cost  ratio  value  which  prevents  the  Energy  Manage¬ 
ment  System  to  adequately  limit  the  batteries  depletion  with  the 
contribution  of  the  Fuel  Cell  power. 

Three  commonly  found  driving  cycles  was  used  to  validate  and 
perform  the  comparative  study.  Compared  with  a  rule-based 
Energy  Management  System,  the  optimal  and  real-time  method 
proposed  in  this  paper  provide  up  to  5%  of  hydrogen  saving  whilst 
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maintaining  the  Battery  Pack  energy  close  to  the  minimum 
prescribed  energy.  The  high-power  demand  reduction  and  the 
depletion  limitation  enable  the  Energy  Management  System  to 
preserve  the  batteries  life. 
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